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ABSTRACT 


A  band  structure  model  is  proposed  in  Fig.  1(d)  which  is 
consistent  with  magneto-Hall  and  magnetoresistance  data  obtained 
for  various  samples  of  HgTe  at  4.2°K.  Using  an  optimization 
routine  it  was  ascertained  that  the  data  could  be  fitted  on  the 
basis  of  a  three-carrier  model  involving  two  sets  of  electrons 
and  one  set  of  holes.  Under  special  crystal  growth  conditions 
the  electron  Hall  mobility  at  4.2°K  of  HgTe  exceeded  640,000 
cm^/v-sec.  •* 
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Hall  Coefficient  and  Transverse  Magnetoresistance 
in  HgTe  at  4.2°K  and  77°K 


Introduction 

While  much  work  has  been  published  in  the  last  few  years 

pertaining  to  the  thermoelectric,  thermomagnetic,  and  optical 

properties  of  HgTe,  there  is  as  yet  no  unanimity  of  opinion 

concerning  the  details  of  the  band  structure  of  this  material. 

In  1961,  it  was  concluded ^  that  HgSe  and  HgSe^  5TeQ  5  are 

sendmetals  and  it  was  suggested  that  HgTe  might  be  likewise. 

This  was  later  corroborated  in  measurements  of  the  Hall  coeffi- 

(2) 

cient  and  of  electrical  resistivity  on  HgTe  and  on  its  alloys, 

at  low  concentration,  with  CdTe.  These  and  related  experiments 

(3) 

were  ultimately  interpreted  in  terms  of  a  band  structure 

(4) 

model  derived  from  that  proposed  originally  by  Groves  and  Paul 
for  grey  tin.  As  illustrated  in  Pig.  1(a),  the  band  structure 
normally  encountered  in  Group  IV  elements  involves  a  conduction 
band  separated  by  a  set  of  Tg  bands  that  are  degenerate  at  the 
origin.  The  dispersion  relations  for  grey  tin  are  shown  in 
Fig.  1(b)  in  terms  of  the  Tg  bands  lying  above  the  Tg  band; 
relative  to  the  diamond  structure,  the  order  of  the  bands  is  now 
inverted . 

The  success  of  the  band  inversion  procedure  in  explaining 
the  semimetallic  properties  of  grey  tin  suggests  that  a  similar 
scheme  might  be  tried  for  HgTe.  The  parent  band  structure  in 
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(b)  GREY  Sn 


Fig.  1.  Various  band  structure  models. 
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this  case,  which  corresponds  to  the  zinc  blende  crystal 

structure,  is  shown  in  Fig.  1(c)  and  is  labelled  as  the  InSb 

structure.  Owing  to  lack  of  inversion  symmetry,  the  energy 

dispersion  relations  now  contain  a  term  linear  to  the  wave 

number  vector,  As  a  consequence,  a  small  split  is  predicted, 

away  from  -£  =  0 ,  in  the  Tg  band  of  Fig.  1(c),  with  shallow 

extrema  located  away  from  the  origin  of  the  Brillouin  zone. 

Apparently,  this  deviation  from  (la)  is  so  slight  that  it  has 

not  been  observed,  even  in  InSb, experimentally .  The  problem  of 

(5 ) 

inverting  the  InSb  model  has  been  considered  by  Gorzkowski  , 
who  indicates  that  a  quantitative  specification  of  e (£)  would 
be  a  complex  undertaking.  In  the  absence  of  a  simple  quantitative 
theoretical  inversion  scheme,  we  therefore  propose  a  band 
structure  model  for  HgTe  ‘in  Fig.  1(d)  that  is  consistent  with 
our  experimental  results . 

If  Fig.  1(b)  represents  the  band  structure  of  HgTe,  then  it 
should  be  possible  to  explain  all  electrical  measurements,  down 
to  liquid  helium  temperatures,  in  terms  of  one  set  of  holes  and 
one  set  of  electrons,  so  long  as  the  Fermi  level  remains  close 
to  the  degenerate  band  edge  energy.  By  contrast,  the  model  re¬ 
presented  by  Fig.  1(d)  predicts  that  two  sets  of  electrons  and  one 
set  of  holes  participate  in  conduction  processes  under  these 

conditions.  This  point  is  brought  out  in  Appendix  A. 

The  interpretation  of  experimental  results  with  respect 

(6 ) 

to  this  problem  is  still  rather  confused:  In  a  review'  of 
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work  completed  through  1964,  it  was  pointed  out  that  a  two- 

carrier  model  is  inadequate  to  account  for  the  transverse 

magnetoresistance  and  magneto-Hall  data  measurements  at  77°K 

and  4 . 2°K  for  a  variety  of  doping  levels,  and  this  was  attributed 

to  an  unusually  complex  valence  band  structure .  More  recent 

work  has  furnished  further  corroboration:  Ivanov-Omskii  and 
(7) 

co-workers  have  reported  systematic  measurements  on  a  series 

of  samples  which  varied  over  a  wide  range  of  net  donor  and 

acceptor  concentration.  Their  attempts  to  fit  the  data  to  a 

(8 ) 

two-carrier  model,  using  the  method  of  Hilsum  and  Barrie  '  ' 

(9) 

proved  unsuccessful.  By  contrast,  Stradling  and  Antcliffe 
have  reported  quantum  oscillations  in  magnetoresistance  at  4.2°K 
in  high  magnetic  fields  on  a  sample  of  HgTe  in  which  the  Hall 
coefficient  reversed  sign  with  increasing  magnetic  field.  They 
interpreted  their  results  in  terms  of  a  two-carrier  model  and 
they  report  an  overlap  energy  in  excess  of  0.0007  eV  between  the 
valence  and  conduction  bands.  No  attempt  was  made  to  determine 

9 

if  their  data  were  also  consistent  with  a  three-carrier  model 
at  low  temperatures.  Delves  investigated  the  properties  of 

single  crystals  of  Hgn  Mn  Te?  from  measurements  at  77°K  for 
small  x  he  provisionally  identified  a  second  set  of  electron- like 
carriers.  On  the  other  hand,  Tovstyuk  and  co-workers ^ 
concluded  that  their  measurements  indicated  the  presence  of  one 
set  of  electrons  and  two  sets  of  holes  in  their  samples. 
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.  (12) 

Kolosov  and  Sharavskii  reach  a  similar  conclusion  on  the  basis 

of  Seebeck  coefficient  measurements  at  188°K  in  a  transverse 
magnetic  field.  Veri4>  bases  the  interpretation  of  his 

electrical  measurements  in  the  range  77  -  373°K  on  a  two-band 
model,  though  he  does  sketch  an  overlap  band  model  similar  to 
that  in  Fig.  1(d).  Finally,  Yamamoto  and  Fukuroi^^)  on  the 
basis  of  the  Shubnikov-de  Haas  effect  conclude  that  HgTe  is 
characterized  by  two  distinct  conduction  bands  involving  electrons 
of  different  mass . 

Part  of  the  difficulty  in  arriving  at  a  clear  decision 
concerning  the  applicability  of  the  various  band  structures  stems 
from  the  fact  that  a  three-carrier  model  yields  very  complex 
expressions  for  the  transport  coefficients^'’^,  so  that  it  be¬ 
comes  difficult  to  make  comparisons  between  theoretical  predictions 
and  experimental  observations.  The  present  work  was  undertaken 
to  clarify  aspects  concerning  the  applicability  of  the  three- 
carrier  model  to  HgTe.  Towards  this  end  we  accumulated  a  large 
body  of  experimental  results  on  the  Hall  coefficient  and  the 
magnetoresistance  of  undoped  HgTe.  Selected  representative  runs 
were  then  subjected  to  optimal  fitting  procedures  by  a  computer 
to  ascertain  whether  a  three-carrier  model  was  indeed  required  to 
fit  the  data,  and  if  so,  whether  one  could  distinguish  between 
the  alternatives  of  requiring  two  sets  of  electrons  and  one  set 
of  holes  as  against  the  inverse  possibility. 


5 


Preparation  of  Materials 


The  three  samples  reported  on  in  this  investigation  were 
prepared  by  the  technique  discussed  in  Ref.  5  With  the  following 
differences  in  detail:  Sample  4A  was  an  as-grown  sample  cut  from 

the  f irst-to-freeze  part  of  an  ingot  grown  at  a  rate  of  3  mm/hr 
from  a  stoichiometric  100  gm  charge  which  had  been  sealed, 

along  with  2.7  mg  of  Al,  in  a  quartz  tube.  The  temperature  of 
the  coolest  part  of  the  tube  in  contact  with  vapor  during 
growth  was  560°C,  corresponding  to  a  mercury  pressure  of  15  atm. 
Sample  300- 2A  was  cut  from  the  f irst-to-freeze  part  of  an  ingot 
grown  at  a  rate  of  3  mm/hr,  from  a  stoichiometric  400  gm  charge 
containing  10.8  gms  of  excess  Hg.  The  temperature  of  the 

coolest  part  of  the  tube  in  contact  with  vapor  during  growth  was 

590°C.  The  as-grown  sample  was  sealed  in  an  evacuated  quartz 
tube  and  the  tube  placed  in  the  first  zone  at  a  temperature  of 
350°C,  while  the  temperature  of  the  coolest  part  or  second  zone 
was  300°C.  The  sample  was  annealed  for  336  hrs  and  quenched  in 
water.  Sample  50  LB  was  cut  3  mm  from  the  first- to- freeze  part 
of  an  ingot,  which  was  grown  from  a  stoichiometric  400  gm  charge 
containing  15  gms  of  excess  Hg .  The  growth  rate  was  0.1  mm/hr,  and 
the  material  was  quenched  in  water.  The  temperature  of  the 
coolest  part  of  the  tube  in  contact  with  vapor  during  growth  was 

590°C.  Since  only  10  gm  of  excess  Hg  was  required  to  make  up  for 

i 

the  loss  of  mercury  to  the  vapor,  this  ingot  crystallized  from  a 
liquid  composition  that  was  slightly  rich  in  mercury. 
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Experimental  Results 

Representative  data  for  the  Hall  coefficient  in  low  magnetic 
field,  for  the  electrical  resistivity  in  the  absence  of  magnetic 
fields,  for  transverse  magnetoresistance,  and  for  Hall  mobility 
at  300,  77  and  4 . 2°K  are  assembled  in  Table  1.  An  interesting 

feature  of  the  data  is  the  unusually  large  value  of  the  electron 
Hall  mobility  of  unannealed  sample  50  LB  at  4.2°K.  The  value  of 
2 

640,000  cm  /volt- sec  is  significantly  higher  than  the  previously 

2 

higher  value  of  360,000  cm  /volt- sec  reported  in  Ref.  6.  In  fact, 
unannealed  specimens  from  ingot  50  have  exhibited  unusually  well 
defined  quantum  oscillatory  effects  in  both  magnetoresistance  and 
magneto-optical  measurements. 

Analysis  of  Results 

The  magneto-Hall  and  magnetoresistance  measurements  taken 

at  77°K  and  4 . 2°K  were  first  analyzed  according  to  the  method  of 

(8) 

Hilsum  and  Barrie.  f  Plots  were  made  of  the  Hall  coefficients 

vs  magnetoresistance  and  of  the  magneto-Hall  and  magnetoresistance 

vs  the  inverse  square  magnetic  field.  These  should  result  in 

straight  lines  on  the  basis  of  a  two-carrier  model;  they  were, 

in  fact,  found  to  exhibit  serious  departures  from  linearity. 

This  clearly  requiren  a  more  sophisticated  analysis  in  terms  of 

at  least  a  three- carrier  model. 

The  transport  theory  involving  such  a  model  was  first 

(17) 

worked  out  by  Willards'on  and  co-workers  for  the  case  of 
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p-type  Ge .  The  final  results  rested  on  two  assumptions,  namely 
that  classical  statistics  was  applicable  and  that  the  mean  free 
path  of  the  charge  carriers  was  independent  of  energy.  Even  so, 
the  final  results  were  rather  complex,  involving  exponential 
integrals  and  error  functions .  In  HgTe  the  situation  is  further 
complicated  by  the  fact  that  for  semimetals  classical  statistics 
is  inapplicable.  Accordingly,  it  was  decided  to  adopt  a  thermo¬ 
dynamic  formulation for  isotropic  materials  in  which  the 
overall  conductivity  o  and  Hall  coefficient  are  expressed  in 

terms  of  the  one-band  contributions  and  (i  =  1,2,3)  . 

That  these  relations  should  apply  to  HgTe  is  substantiated  by 

unpublished  work  of  the  authors  on  the  angular  dependence  of 

magnetoresistance  of  oriented  single  crystals  of  HgTe  in  high 

magnetic  fields;  no  significant  anisotropy  effects  were  encountered 

even  in  the  mixed  conduction  region . 

For  ease  of  interpretation,  the  substi tut ions  =  n^qu^ 

and  =  1/n^e  were  used  in  the  overall  expressions  for  and  a  . 

Here,  e  is  the  electronic  charge,  n^  the  carrier  density,  and 

the  carrier  mobility  of  the  ith  species.  In  adopting  these 

expressions  we  neglect  the  variation  of  ^  with  the  magnetic 

field,  H,  which  is  usually  considerably  less  than  the  maximum 

value  of  50%.  We  also  neglect  the  dependence  of  a.  on  H,  which 

(18  20) 

is  usually  considerably  less  than  the  maximum  factor  of  2.4.  *  ' 

Since  changes  in  ^  and  a  of  many  orders  of  magnitude  are  observed. 
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it  is  believed  that  the  conclusions  reached  in  this  publication 

are  not  altered  by  use  of  these  simplifications. 

(19) 

As  Fischer  wfis  first  to  show,  in  the  approximation 

scheme  discussed  above,  the  dependence  of  the  overall  resistivity  p 
and  of  the  overall  Hall  coefficient  (\_  on  H  may  be  expressed  as 


P  (x) 


1  A  +  Bx  +  Lx2 

e  ~2 - 2 

A  +  Dx  +  Mx 


(1) 


6L(x) 


1  E  +  Fx  +  Nx2 

e  2  2 

e  A^  +  Dx  +  Mx 


(2) 


where,  with  e  in  coulombs,  x  e  10  H  (H  in  gauss),  and  where 

A,  B,  D,  E,  F,  L,  M,  and  N  are  eight  nonlinear  functions  of 

the  three  charge  carrier  densities  n^  (in  cm  )  and  of  the  three 

2 

mobilities  p^  (in  cm  /volt- sec)  participating  in  conduction. 

The  parameters  n^  and  p^  were  determined  with  the  aid  of 
a  computer  program  described  in  further  detail  in  Appendix  B. 

In  essence,  the  experimentally  determined  p(0)  and  (H^)  were 
used  as  input  variables  for  a  program  which  yielded  values  of 
n^  and  p^  that  provided  an  optimal  fit  of  Eq.(2)  to  the  experi¬ 
mental  Hall  measurements.  With  the  n^  and  p^  determined  in  this 
manner,  the  magnetoresistance  was  then  calculated  from  Eq. (1) 
according  to  the  relation  p(H)/p(0)-l. 

The  parameters  determined  in  this  fashion  for  three  repre¬ 
sentative  specimens  are  listed  in  Table  2.  The  following  points 
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Fig.  2.  Theoretical  and  experimental  Hall  coefficient  of  HgTe  at  4.2°K. 
Samples  50  LB,  4-A  and  300-2A. 
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R  (cm3/c) 


should  be  noted:  (1)  The  signs  for  the  n^  and  were  determined 
by  the  computer;  from  the  tabulation  one  concludes  that  an 
optimal  fit  of  the  data  requires  that  two  of  the  three  carriers 
be  electron-like .  (2)  A  variation  of  any  of  the  carrier  density 
cr  mobilities  by  ±10%  about  the  listed  value  resulted  in  a 
noticeably  poorer  fit  of  the  theory  to  the  data,  as  judged  from 
least-squares  calculations.  This  was  also  the  case  when  one  or 
more  of  the  parameters  were  simultaneously  altered  by  ±10%.  The 
values  listed  in  Table  2  as  providing  the  optimal  fit  are  thus 
known  with  reasonable  precision. 

A  comparison  between  the  calculated  and  the  observed  magneto- 
Hall  effects  is  provided  in  Fig.  2.  Correlations  between  the 
results  of  Table  2  and  the  observations  in  Fig.  2  are  readily 
apparent:  Sample  50  LB  exhibits  a  negative  Hall  coefficient 

because  the  two  sets  of  electrons  outweigh  the  holes  in 
concentration.  In  passing  one  should  note  the  distinct  maximum 
(actually,  a  minimum)  in  the  magneto-Hall  curve  for  this  sample; 
on  the  basis  of  elementary  arguments  it  may  be  shown  that  this 
cannot  occur  in  a  two-carrier  model  under  the  assumption  that 
the  one  carrier  6^  and  are  independent  of  H.  In  samples  4A 
and  300-2A  the  sign  reversal  of  tK.  in  passing  from  low  to  high 
magnetic  fields  reflects  the  preponderance  in  concentration  of 
low-mobility  holes  over  intermediate  and  high-mobility  electrons. 

As  is  evident  from  Fig.  3,  there  is  considerable  discrepancy 
between  the  calculated  and  the  observed  magnetoresistance  values. 
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Fig.  3.  Theoretical  and  experimental  magnetoresistance  of  HgTe  at  4.2°K.  Samples  50  LB,  4 -A 
and  300-2A. 


although  the  general  trend  of  the  experimental  results  is 

correctly  brought  out  in  the  calculations.  The  fit  could  be 

considerably  improved  by  using  both^(H^)  and  p(H^)  as  input 

parameters  to  the  computer  program.  The  discrepancies  are  due  in 

part  to  the  simplifying  assumption  which  neglects  the  variation 

of  ^  and  on  H.  It  should  be  noted  that  in  sample  50  LB  the 

magnetoresistance  attained  values  in  excess  of  600  without 

showing  any  saturation  effects  at  110  kG. 

The  results  cited  above  are  consistent  with  the  band  model 

sketched  in  Fig.  1(d).  As  has  been  pointed  out  in  Ref.  6,  it 

is  generally  rather  difficult  to  determine  the  thermal  overlap 

energy  with  precision.  A  lower  limit  to  this  quantity  is 

calculated  from  the  value  of  n^  for  sample  50  LB?  here, 

where  E„  is  the  Fermi  level  relative  to  the  bottom  of  the 
F 

conduction  band.  Using  an  effective  electron  mass  of  m/mQ  =0.023, 

( 16} 

as  determined  from  the  direct  gap'  '  energy  value  of  E^  =  0.283  eV 
and  from  the  Kane  parameter  (p  =  8.3  x  10~®  eV-cra,  one  finds  that 


0.0028  eV? 


(3) 


an  upper  limit ^  to  E^  of  0.02  eV  is  set  by  calculations  involving 
the  two-band  model . 
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Conclusions 


Magneto-Hall  and  transverse  magnetoresistance  data  on  a 
large  number  of  HgTe  single  crystals  indicate  that  a  three-carrier 
model  involving  two  sets  of  electrons  and  one  set  of  holes  should 
be  used  in  the  interpretation  of  the  data.  These  findings  are 
consistent  with  an  inverted  band  structure  model  shown  in  Fig.  1(d); 
the  lower  limit  to  the  overlap  energy  E^.  is  0.0028  eV.  An  unusually 
high  electron  mobility  of  640,000  cm^/V-sec  was  obtained  at  4 . 2°K 
as  a  result  of  special  conditions  of  crystal  growth. 
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Appendix  A  -  The  Quartic  Band 


In  this  appendix  we  describe  briefly  the  electrical  charac¬ 
teristics  to  be  associated  with  a  quartic  band  model  of  the 
form 


e 

kT 


x 


(A-l) 


which  is  sketched  in  Fig.  4.  Here  e  is  the  energy,  k  is 
Boltzmann's  constant,  A  and  B  are  constants,  and  is  the 

radial  component  of  the  wave  number  vector.  It  is  seen  that  a 
band  of  the  form  (A-l)  simulates  one  of  the  ro  bands  in  the 

O 

proposed  band  scheme  for  HgTe,  Fig.  1(d). 

It  is  important  to  note  that  for  such  a  band  with  spherical 
symmetry  the  current  density  may  be  split  into  two  contributions, 
namely 


J  = 


+ 


Zev 

'V 


f,  -fe  ^  sinO  d0  d0. 


+  J  q 


(A- 2) 


where  e  is  the  electronic  charge,  Z 


±1,  v  =  IT1  e,  f-L  the 


deviation  of  the  nonequilibrium  distribution  function  from  the 
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BAND  STRUCTURE  CHARACTERIZED 

€  4  2 

BY  —  =  x  =  -Ak  +  Bk 
kT 

Fig.  4.  The  quartic  band  structure  model. 
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Fermi-Dirac  distribution  function  f  =  [exp  (x-t)) +1]_\  t)  the 
reduced  Fermi  level  relative  to  the  appropriate  band  edge;  6 

and  0  are  the  appropriate  angular  coordinates  in  _  space,  and 

h.  ~ 

quantities  v  and  are  functions  of  x  alone  and  are  thus  inde¬ 
pendent  of  9^  ,  0^  .  One  may  therefore  treat  the  band  in  Fig.  4 
as  a  combination  of  two  sub-'bands,  each  of  which  may  be  considered 
separately . 

Focusing  on  the  inner  band,  0  <  ,  and  switching  to 

the  formulation  of  transport  effects  in  terms  of  irreversible 
thermodynamics  we  now  write: 
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where  V  =  v  (£/e),  £  being  the  Fermi  level  relative  to  the 
A  A 

electron  at  xest  in  free  space,  and  where  a  (-H  )  =  g  (H  ), 

yx  z  xy  z 

2  2 

axx  (Hz  )  =  °yy  )  are  the  components  of  the  conductivity  tensor. 

Setting  J^y  =  0  and  eliminating  V  from  (A- 3)  we  then  obtain 
(V  T  =  0,  J,y  =  0  ) 


ayx^  axyayx 


axxayy  ^ 


(A- 5) 


As  the  band  under  consideration  is  filled  with  charge  carriers. 
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the  denominator  in 


thereby  increasing  £  from  0  towards  ^5.  , 

(A- 5)  remains  unaltered  in  sign.  The  sign  of  is  thus  deter - 

(20) 


mined  by  that  of  a 


yx 


As  shown  elsewhere 


for  bands  of 


spherical  symmetry,  and  for  cases  where  the  Boltzmann  transport 
equation  in  the  relaxation  time  formalism  applies 


a  =  e 
yx 


rK  {*/2  r2lx  vy 

J  J  J  1  +  oi2r2 

°  "  7T/2  ° 


T2  V  sin0k  d  *r  d0*  dt' 

de  ' 


(A- 6) 


where  r  is  the  relaxation  time,  v’  =  Vv  e>  and  where 

y  Ky 

co  =  -eHz/m*c.  The  sign  of  <jyX  is  thus  governed  by  that  of  the 
first  der ivative  effective  mass 

1/m*  =  (l/4r)  (de/dfc  r).  (A-7) 

With  being  the  radius  vector  relative  to  the  origin  of  coor¬ 

dinates  at  the  point  of  spherical  symmetryj  the  effective  mass, 
according  to  (A-l),  remains  positive  for  0  <  A  r  <  ^  o’  but  is 
negative  for  Jk.  >  J .  It  follows  that  insofar  as  the  Hall 
coefficient  is  concerned,  the  inner  portion  of  the  band  in  Fig.  4 
is  always  populated  with  electron- like  carriers,  regardless  of  the 
degree  of  filling,  while  the  outer  portion  of  the  band  in  Fig.  4 
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is  populated  with  holes.  It  is  of  interest  that  a  similar  argu¬ 
ment  pertaining  to  the  Seebeck  coefficient  leads  to  the  conclusion 
that  the  thermoelectric  power  changes  sign  as  the  inner  band  is 
progressively  filled  with  carriers. 

These  findings  are  substantiated  by  a  detailed  study  in 
which  (A-l)  was  solved  for  £^(x)  and  the  results  used  in  the  ex¬ 
plicit  evaluation  of  transport  integrals  according  to  the  pro¬ 
cedures  of  Ref.  20. 
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Appendix  B  -  Optimization  Routine 


The  charge  carrier  densities  and  mobilities  listed  in  Table 

2  were  determined  by  an  optimization  routine  due  to  Nelder  and 

(21)  .  ... 

Mead.  In  essence,  the  program  is  designed  to  minimize  the 

function 
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(B-l) 


where  the  subscript  E  and  T  designate  respectively  the  Hall  coef¬ 
ficients  determined  experimentally  and  those  calculated  according 
to  Eq. (2) .  In  Eq. (B-l)  the  variable  n^  has  been  eliminated  by 
use  of  Eq.(l),  rewritten  as 


A  =  n.j y1  +  n 2/i2  +  n3jLt3  =  1/ep  (o) 


(B-2 ) 


The  method  requires  an  initial  guess  of  n n^,  y^,  y  ,  y^ 
which  locates  a  point  in  fdve-space;  the  computer  then  generates 
five  additional  point's  in  the  vicinity,  according  to  a  pre¬ 
scribed  routine,  and  determines  the  center  of  gravity  of  the  re¬ 
sulting  convex  polyhedron  in  five-space,  called  a  simplex.  The 
function  0  is  then  calculated  according  to  Eq. (B-l)  for  each  of 
the  vertices  of  the  simplex;  the  vertex  (extremum)  corresponding 
to  the  highest  0  is  singled  out,  and  a  new  point  is  located  at 
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a  prescribed  distance  along  the  line  passing  through  the  vertex 
in  question  and  the  center  of  gravity.  The  old  extremum  is  now 
discarded,  and  a  new  simplex  generated  from  the  new  point  and 
the  remaining  vertices  of  the  original  simplex.  The  procedure 
is  repeated  until  the  final  simplex  consists  of  six  points 
spaced  very  closely  together.  The  center  of  gravity  now  cor¬ 
responds  to  the  optimal  values  of  the  parameters;  however,  this 
set  under  certain  conditions  may  correspond  to  a  local  minimum 
in  0  rather  than  to  the  true  minimum.  To  guard  against  this 
possibility,  each  of  the  points  of  the  final  simplex  is  distended 
in  turn  by  5  to  10%  from  its  value,  and  the  process  of  finding 
a  minimum  in  0  is  repeated. 

The  advantage  in  the  approach  lies  in  the  fact  that  it  is 
readily  programmed  in  Fortran,  that  it  is  free  from  erratic  loops, 
and  that  it  does  not  require  programming  of  the  derivatives  of 
0 .  The  total  machine  time  required  to  run  through  a  set  of  ex¬ 
perimental  points  on  a  SDS  930  computer  was  of  the  order  of 
five  minutes . 
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